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ABSTRACT 

In a universe with inhomogeneous reionization, the ionized patches create a second 
order signal in the cosmic microwave background polarization anisotropy. This signal 
originates in the coupling of the free electron fluctuation to the quadruple moment of 
the temperature anisotropy. We examine the contribution from a simple inhomogeneous 
reionization model and find that the signal from such a process is below the detectable 
limits of the Planck Surveyor mission. However the signal is above the fundamental 
uncertainty limit from cosmic variance, so that a future detection with a high accuracy 
experiment on sub-arcminute scales is possible. 

Subject headings: cosmic microwave background - cosmology:theory 

1. Introduction 

The Microwave Anisotropy Probe and the Planck Surveyor (MAP; Planck Surveyor; Bersanelli 
et al. 1996) will provide a precise measurement of the cosmic microwave background (CMB) tem- 
perature anisotropy. While existing measurements of the polarization of the CMB radiation only 
give crude upper limits to its anisotropy (Wollack et al. 1993; Partridge et al. 1997), forthcoming 
experiments will be more sensitive to the power in the polarization (Staggs, Gundersen & Church 
1999). MAP is expected to make only a statistical detection of the polarization anisotropy, while 
Planck will measure the polarization power spectrum to high accuracy with an average pixel sensi- 
tivity for the polarization fluctuation of around AT/T = 5 x 10~^ (Planck Surveyor). Since Planck 
will measure the power spectrum up to arcminute scales (Planck Surveyor), second order effects in 
the fluctuation of the CMB radiation could become important. The study of second order effects is 
well established for the temperature anisotropics (Sunyaev & Zel'dovich 1970, 1980; Kaiser 1984; 
Ostriker &: Vishniac 1986; Vishniac 1987; Efstathiou 1988; Dodelson & Jubas 1995; Aghanim et 
al. 1996; Gruzinov & Hu 1998; Peebles & Juskiewicz 1998; Knox, Scoccimarro & Dodelson 1998; 
Haiman & Knox 1999), however for the polarization power spectrum these effects have not been 
examined extensively in the past (Sunyaev & Zel'dovich 1980; Efstathiou 1988) and are only the 
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subject of very recent investigations (Seshradi & Subramanian 1998; Hu 1999). It is well known 
that the universe is ionized at least out to redshifts of z ~ 5 (Gunn & Peterson 1965) and any 
realistic model of how this reionization might have taken place is thought to be inhomogcneous. 
The hot ionized gas interacts with the CMB photons and results in features in the anisotropy power 
spectrum. In this Letter we will focus on the second order effects from inhomogcneous reionization. 
The first order effect of reionization is an enhanced polarization anisotropy on large scales and 
a suppression on small scales (Zaldarriaga 1997). The second order effect due to homogeneous 
reionization is a Ostriker-Vishniac-type effect for the polarization, where the fluctuation in the free 
electron density follows the linear density variations of the overall matter (Efstathiou 1988; Seshradi 
&: Subramanian 1998). The second order contribution we are going to discuss in this Letter is very 
similar to this effect, only that the source of the fluctuation in the free electron density is different. 



2. Second order polarization anisotropy from inhomogeneous reionization 

The source of reionization is thought to be the UV radiation of early objects like quasars 
and proto-galaxies as hosts of an early generation of stars. The nuclear and gravitational energy of 
these objects is transformed into radiation which subsequently ionizes the hydrogen in spheres which 
surround them (Tcgmark, Silk &: Blanchard 1994; Roes 1996; Aghanim et al. 1996; Haiman & Loeb 
1997. 1998; Loeb 1997; Silk k Rccs 1998; Haiman & Knox 1999). One way to study the consequences 
of inhomogcneous reionization is to use effective models which describe the distribution of ionized 
regions by a small number of free parameters (Gruzinov & Hu 1998; Knox, Scoccimarro & Dodelson 
1998; Haiman & Knox 1999). We adopt the model by Gruzinov & Hu (1998), which describes 
inhomogcneous reionization as a set of uncorrelated patches of a certain fixed size R. The number 
density of these patches grows with time and finally the whole universe becomes reionized in a 
homogeneous way. More realistic models include varying patch sizes (Aghanim et al. 1996) and 
contain correlations of the ionized regions. It turns out that these correlations lead to a somewhat 
different signal on smallest scales for the temperature anisotropies (Knox, Scoccimarro & Dodelson 
1998). However, the naive uncorrelated model by Gruzinov &; Hu (1998) gives a good estimate of 
the effect on the CMB (Haiman & Knox 1999). 

We define the ionization fraction Xc to be the ratio of the number density of free electrons Uq and 
the overall (free and bound) number density of electrons n, i.e. = Ue/n. Since we want to study 
the effects of inhomogeneous reionization on the CMB to second order, we not only need the mean 
background ionization fraction Xe, we also have to know the variance of the distribution. Gruzinov 
Sz Hu (1998) give the second moment to be 

(^l-X2)^ 

(a;e(?7l,Xi)a;e(?72,X2)) = a;e(?7min)a;e(??max) + a;e(??min) [1 - a;e(?7max)] e 2Jl^ , (1) 

where i] = J dt/a is the conformal time, with a the scale factor, rjram = 111111(171,772) and 77max = 
max(?7i, 772). The correlation (1) drops off exponentially if the distance between two points is 
larger than the size of a patch R. If 771 or 772 is in one of the homogeneous regimes (nearly 
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neutral or complete ionization) the correlation is just the product of the mean values. In the 

following we will describe the inhomogeneity as the fluctuation of the free electron number density 
(5e(??,x) = (ne(??,x) — nc{r])) / nc{rf) . We study a universe with a 5% baryon content, a Hubble 
constant of Hq = 100/ikm/sec/Mpc, with h = 0.5, critical matter density and no cosmological 
constant, but the results are easy to generalize to an open or A universe. The matter fluctuations 
are taken to be adiabatic with an initial spectral index of n = 1.0. 

We are not including tensor perturbations and therefore the magnetic component of the polarization 
fluctuation is zero (Zaldarriaga & Seljak 1997). It is then sufficient to study the perturbations in 
the Stokes parameter Q. If we assume the wave is traveling into the z direction with and 
Ey the amplitudes of the electric field in x and y direction respectively, this parameter is given 
by Q = {Ex) — i^y)- The first order Boltzmann equation for the fluctuations in Q is (Bond & 
Efstathiou 1987; Zaldarriaga & Seljak 1997) 

Aq + jidiAq = n^aTa (^-Aq + ^ [1 - P2{fi)] , (2) 

with 7j the direction of the photon momentum, -P2(m) tke second Legendre polynomial, 11 = 
At2 + Aq2 + Aqo the polarization tensor and the overdot refers to the derivative with respect to 
conformal time. We expand in this equation the free electron density to first order fluctuations 
ne(?7,x) = ne(r/) [1 + (5e(r/,x)] and collect all second order terms which involve the free electron 
fluctuation Sc. Since we are only interested in the consequences of inhomogeneous reionization 
we do not include other second order contributions. The integral solution of this second order 
contribution in Fourier space is then 

Ag)(%,k,/x) = -(1 - /.2) gir,o,v)e"'^'^-'°^^S{v,\^)dv , (3) 

where /i = cos^? = jiki/k and rjo is the conformal time today. We have neglected in this solution 
the couplings of the first order polarization fluctuations to 6^, since the first order temperature 
quadrupole will dominate these terms (Seshradi &; Subramanian 1998). The homogeneous back- 
ground ionization history is encoded in the visibility function g{rjQ,r)) = rexp{— r(?7o) +t{ii)} with 
the differential optical depth f = x^naT^a, where cjt is the Thomson scattering cross section. The 
source S of the second order fluctuation is the mode coupling term between the fluctuation in the 
free electron density and the first order quadrupole fluctuation of the temperature, 

S{V,^) = J 5,{v,k-p)Al^l{r,,p)d^p^j^Sc{v,k) J Ag(r?,p)d=^p. (4) 

The approximation in (4) needs explanation and we follow the argument of Seshradi & Subramanian 
(1998). The free streaming solution of the flrst order quadrupole of the temperature anisotropy is 
proportional to the second spherical Bessel function, A^2('?) P) c< jiipiv — ''?rec])) with r/rec the time 
of recombination (Hu &; Sugiyama 1995). The second spherical Bessel function can be approximated 
by a Gaussian with a peak around p ^ = 3.345/(7/ — r/rec) and since we are only interested in 
reionization times below z = 100 we get po < (250/i""^Mpc)""^. The typical size of a reionized patch 
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is of the order 10/i~^Mpc, which defines the scale where the free electron fluctuation 6e varies the 
most. Therefore 6c is nearly constant where the quadrupole has the largest contribution to the 
integral in (4), i.e. (5c(??, k — p) ~ 5c(??, k — po) in the relevant integration range. Furthermore the 
dominant /c-range is of the order (10/i~^Mpc)~^ and much larger than pQ, therefore we can write 

5e(77,k-po) « (5e(77,k). 

To calculate the power spectrum we expand the fluctuation in Q into spin-2 spherical harmonics 
and calculate the correlator of the expansion coefficients (Zaldarriaga & Seljak 1997). The Stokes 
parameter Q is not invariant under rotation and therefore dependent on the coordinate system we 
choose. However we can construct an invariant quantity by applying a spin raising operator on Q 
(Zaldarriaga & Seljak 1997). In our case the resulting quantity is just the electrical field E-type 
component of the polarization fluctuation. To calculate the two-point function of the expansion 
coefficients we need to know the correlator of the source S{r],]s.), which is 

oo 

(<S(r?i,ki)<S(r?2,k2)) «i 47r {Se{viM)Se{V2M)) J dpp^Pi{p)A!}l{ni,p)A!}l{n2,p) , (5) 



where we have applied the approximation from equation (4). Further we exploit the fact that 
the first order temperature quadrupole anisotropy is uncorrelated to the free electron fluctuation 
from inhomogeneous reionization. This is not necessarily true in a universe where the patches are 
correlated but the expression can be worked out as long as the underlying fluctuations are Gaussian. 
The quadrupole fluctuation iv^ P) equation (4) is a random variable with amplitude and phase 
depending on the initial perturbations V'i(p) with A^2('7) P) = V'i(p)^T2('?'^') (V'i(Pi)V'i(P2)) = 
Pi{p)S^^^ (Pi — P2) (Ma &; Bertschinger 1995). The correlator in the free electron fluctuation is given 
by (1) and the integral in (5) is the correlator of the first order quadrupole fluctuation at unequal 
times. The factor Pi{p) is the COBE normalized initial power spectrum in the metric fluctuations 
and is given in inflationary, adiabatic models by a power law Pi{p) oc jo"^^, with n the spectral 
index. The second order anisotropy power spectrum for the E-mode polarization is then given by 

<^e} = i(27r)^/^i?^§^ J dkdriidri2g{vo,rii)9irio,m)Qp{vi,V2)Iim,m)x 



k"^ exp 



with 1(771,772) = Xg^{riiain) — 1, jlix) the spherical Bessel functions, Xi = k{r]o — rji) and the 
correlation in the quadrupole at unequal times Qp{rii,r]2) = / dpp'^Pi{p)A^l{r)i,p)A^l{r]2,p), i.e. 
Qp{r]o,i]o) oc Ct,2- We should notice that the limits of the time integrals are the start and end of 
the reionization process, i.e. we integrate over the time during which the inhomogeneities appear. 
For the background reionization history we assume that the mean ionization fraction x^ grows from 
zero to unity between the redshifts z* — 6z* /2 and z* + 6z* /2, which we will choose appropriately. 
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3. 



Results 



We have calculated (6) with a modified version of CMBFAST (Seljak & Zaldarriaga 1996). The 
COBE normalized polarization power spectrum is shown in figs.l and 2 for different parameters. 
In fig.l the background reionization is given by z* = 50 and 5z* = 20. The dotted line refers to 
the first order contribution. Wc recognize the feature on large scales which appears because of the 
homogeneous reionization background. The thick long-dashed, solid and short-dashed lines refer 
to the second order contribution due to inhomogeneous reionization. The long-dashed line is a 
model with a patch size of i? = 20/i~^Mpc, the solid line for R = 10/i^^Mpc and the short-dashed 
line for R = 5^~^Mpc. 

To understand the behavior of the second order effect we will perform two approximations to 
the integral in (6). First we realize that the expression g{r]Q,rji)g{r]Q,r]2)I{rji,rj2), for reasonable 
parameters of the reionization history, is a 2-dimensional function with a narrow peak. This allows 
us to perform the two time integrations in (6) and calculate the integrand at a certain time rj* 
multiplied by an area factor {Srj)'^. The second feature we exploit is that the spherical Bessel 
function ji{x) has a tight peak at I = x for large multipole moments I, so we can approximate it 
with a (5-function. Therefore we get 



with 9o = R/{r]o — r]*) as the angular size of a patch as seen today. The only unknown quantity 
in this expression is the effective time period Srj. For the case with R = lO/i^^Mpc we get a 
good fit with 5r] = 6.5/i~^Mpc, which is in the range of the width of the narrow peaked function 
^1)5(^0) ??2)-^(??i, ??2)- We recovered a similar behaviour for the whole range of reasonable 
reionization parameters. Although we can not predict the amplitude of the second order power 
spectrum in (6) analytically, the shape is well approximated by /(/ -|- l)C^] oc /^e~'^®o/2j^(|). This 
describes essentially an rise with some cut-off around the scale I « \/2/0o and we see in 

fig.l how the peak of the signal moves to larger multipoles I when we decrease the patch size R. 
We also recognize that the amplitude is proportional to 9q and therefore the larger the patch size 
R the larger the power in the second order anisotropy from inhomogeneous reionization. 
In fig. 2 we have plotted the results for a model with a more realistic background reionization history 
(Haiman &; Knox 1999). In this case reionization takes place around z* = 10 with a time period 
of Sz* = 5. The long-dashed line corresponds to a patch size oi R = 10/i~^Mpc, the solid line to 
R = 5/i~^Mpc and the short dashed-line to R = l^~^Mpc. One clearly sees that the power in 
the anisotropy for this late time reionization is much smaller than in the case for earlier times in 
fig.l. This is because the visibility girjo, r/*) is much smaller for the short and late time reionization 
phase. Again the the peak moves to the right when the patch size is decreasing. For inhomogeneous 
reionization with correlated patches we expect the same behavior as for the temperature anisotropy 
power spectrum as discovered by Knox, Scoccimarro Sz Dodelson (1998). In these models the power 
in the anisotropies is much wider distributed over the multipole moments I than for an uncorrelated 
scenario, but the magnitude is very similar. 



(Z-2)! 



elr^Qpir]*,r]*)I{ri*,r]*)g\rio,r]*){Sr])''e-^jl{l) 



(7) 
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4. Conclusion 

In figs.l and 2 we find that the second order signal dominates over the first order signal 
only for very large multipoles. Even for an unrealistic reionization model like in fig.l the second 
order contribution is relevant only on scales smaller than 5 arcsec. One might hope, that once 
the cosmological parameters arc estimated by the first order temperature power spectra and other 
experiments, one can reveal the nature of the second order effects. The Ostriker-Vishniac effect 
for polarization is of the same order or smaller than the signal from inhomogcneous reionization, 
dependent on the ionization parameters (Seshradi & Subramanian 1998). However this effect is 
completely determined by the linear power spectrum and therefore can be removed like the first 
order contribution. 

At the present polarization has not been detected in the CMB and the measurements give only 
crude upper limits (Staggs, Gundersen & Church 1999). We have given the 95% confidence upper 
limits from the Saskatoon anisotropy experiment (Wollack et al. 1993) and the VLA 8.4 GHz CBR 
project (Partridge et al. 1997) as a circled and a diamond point in figs.l and 2. The most accurate 
future experiment which will measure the polarization anisotropy is the Planck Surveyor. Its High 
Frequency Instrument (HFI) is expected to measure the polarization to high accuracy in its 143 
and 217 GHz channels. The average sensitivity AT/T to the linear polarization per pixel and the 
angular resolution is 3.7 x 10~^ and 8.0 arcmin, and 8.9 x lO^*' and 5.5 arcmin, respectively (Planck 
Surveyor). We have calculated the expected polarization signal sensitivity due to cosmic variance, 
beam size and instrument noise with the methods described in (Knox 1995; Bond, Efstathiou & 
Tegmark 1997; Bond, JafFe & Knox 1998). In figs.l and 2 we show the sensitivity histograms, 
resulting from a logarithmic bin by weighted average. One recognizes that the noise levels are 
all above the second order signal from inhomogcneous reionization, so that even with a fairly 
wide binning strategy, the sensitivity of the Planck Surveyor is not large enough to reveal these 
signals. The dashed horizontal lines in the histograms in figs.l and 2 are the logarithmically binned 
uncertainty contributions just from cosmic variance, which are given by (AC;)^ = 2Cf /{2l + 1), 
with full sky coverage. Cosmic variance is the fundamental uncertainty limit and describes the fact 
that we can only observe one universe with only 21 + 1 modes of a certain multipole moment I. 
It is clear from figs.l and 2 that the second order contribution from inhomogcneous reionization 
is above these limits for large multipoles. Therefore a high accuracy polarization measurement on 
sub-arcminute scales could reveal such a signal. But this depends on how well one can remove 
polarization foregrounds and the magnitude of other foreground-type second order contributions, 
like the Sunyaev-Zel'dovich effect (Hu 1999). However, there might be the possibility to disentangle 
all these effects by exploiting small scale measurements of the matter distribution and the CMB 
anisotropics, including the information from polarization. 

The author gratefully acknowledges helpful conversations with N. Aghanim, A. Albrecht, 
L. Knox, A. Lewin, J. Magueijo and M. Zaldarriaga. Further I thank U. Seljak and M. Zal- 
darriaga for the use of CMBFAST. JW is supported by the German Academic Exchange Service 
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Fig. 1. — The CMB polarization anisotropy power spectrum for a model with an effective reioniza- 
tion time z* = 50 and Sz* = 20. The dotted line is the first order contribution. The second order 
signal from inhomogeneous reionization is the long-dashed line for a patch size of i? = 20/i-iMpc, 
the solid line for R = 10/i~^Mpc and the short-dashed line for R = 5/i^^Mpc. The diamond and 
circled data points are the 95% confidence upper limit from the Saskatoon anisotropy experiment 
and the VLA 8.4 GHz CBR project, respectively. The histogram shows the logarithmically binned 
uncertainty limit from the polarization measurement with the Planck Surveyor. The horizontal 
dotted lines in the histogram correspond to the uncertainty levels due to cosmic variance only. 
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Fig. 2. — The same plot as figure 1 but for an effective reionization at z* = 10 and 5z* = 5. 
Tfie patcfi sizes are R = 10/i~^Mpc for the long-dashed line, R = 5h~^Mpc for the solid line and 
R = lh~^Mpc for the short-dashed line. 



